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After reviewing the basic relevant properties of stationary stochastic processes (SSP), defining basic 
terms and quantities, we discuss the properties of the so-called Harrison-Zeldovich like spectra. 
These correlations, usually characterized exclusively in fc-space (i.e. in terms of power spectra 
P{k)), are a fundamental feature of all current standard cosmological models. Examining them 
in real space we note their characteristics to be a negative power law tail ^(r) ~ — J'"^, and a 
sub-poissonian normalised variance in spheres u^{R) ~ i?~'*ln_R. We note in particular that this 
latter behaviour is at the limit of the most rapid decay (~ _R~^) of this quantity possible for 
any stochastic distribution (continuous or discrete). This very particular characteristic is usually 
obscured in cosmology by the use of Gaussian spheres. In a simple classification of all SSP into three 
categories, we highlight with the name "super-homogeneous" the properties of the class to which 
models like this, with P(0) = 0, belong. In statistical physics language they are well described as 
glass-like. They do not have either "scale-invariant" features, in the sense of critical phenomena, 
nor fractal properties. We illustrate their properties with some simple examples, in particular that 
of a "shuffled" lattice. 
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I. INTRODUCTION 

In standard theories of structure formation in cosmol- 
ogy the density field in the early Universe is described 
as a perfectly homogeneous and isotropic matter distri- 
bution, with superimposed tiny fluctuations character- 
ized by some particular correlation properties (e.g. [0|). 
These fluctuations are believed to be the initial seeds 
from which, through a complex dynamical evolution, 
galaxies and galaxy structures have emerged. In partic- 
ular the initial fluctuations are taken to have Gaussian 
statistics and a spectrum which is exactly, or very close 
to, the so-called Harrison-Zeldovich (hereafter H-Z) [HJ|] 
or "scale- invariant" power spectrum (hereafter PS). Be- 
cause fluctuations are Gaussian, the knowledge of the PS, 
or its Fourier conjugate, the real space correlation func- 
tion, gives a complete statistical description of the fluctu- 
ations. The H-Z type spectrum was flrst given a special 
importance in cosmology with arguments for its "nat- 
uralness" as an initial condition for fluctuations in the 
framework of the expanding universe cosmology, and it is 
in this context that the use of the term "scale-invariant" 
to designate it can be understood. It subsequently gained 
in importance with the advent of inflationary models in 
the eighties, and the demonstration that such models 
quite generically predict a spectrum of fluctuations of 
this type. Since the early nineties, when the COBE ex- 
periment [Q measured for the first time fluctuations in 
the temperature in the Cosmic Microwave Background 
Radiation (hereafter CMBR) at large scales, and found 
results consistent with the predictions of models with a 
H-Z spectrum at such scales, the H-Z type spectra have 
become a central pillar of standard models of structure 
formation in the Universe. The aim of the present paper 



is twofold. Firstly, to clarify the statistical properties in 
real space of these distributions, which have been almost 
completely overlooked in the literature on the subject. 
And secondly, through this discussion, to relate and com- 
pare this model of the primordial Universe to correlated 
systems encountered in statistical physics. We attempt 
to make the paper as self-contained as possible, and not 
excessively technical in its discussion either of cosmology 
or statistical concepts, in the hope that it may be easily 
accessible to both cosmologists and statistical physicists. 

The H-Z spectrum arises in cosmology through a par- 
ticular condition applied to perturbations of Friedman- 
Robertson- Walker (FRW) models, which describe a ho- 
mogeneous Universe in expansion. This condition - com- 
monly referred to in cosmology as "scale invariance" of 
the perturbations - gives rise to a spectrum (commonly 
called the "scale-invariant" perturbation spectrum) with 
P{k) ~ fc at small k. All current standard cosmological 
models of structure formation in the Universe assume a 
spectrum exactly like this, or close to it, as initial con- 
dition for perturbations in the Universe. In such mod- 
els there is at any time a finite scale corresponding to 
the causal horizon, which increases with time, and below 
which causal physics can act to modify the spectrum. 
This causal physics depends, in general, on the details of 
the model, i.e. on the nature of its content in matter and 
radiation (or other forms of energy) , until a characteristic 
time (the time when matter and radiation have compa- 
rable densities), after which purely gravitational evolu- 
tion takes over. There are many variants on standard 
cosmological models e.g. "Cold Dark Matter" (hereafter 
CDM), "Mixed Dark Matter" (hereafter MDM), or the 
currently favored one with a non zero cosmological con- 
stant (hereafter ACDM) , each of them leading to a differ- 
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ent form for the spectrum at smaller scales (i.e. large k) 
which can be calculated. In CDM models (in which the 
predominant massive component driving collapse under 
gravity is Cold Dark Matter, "cold" in the sense that the 
particles have little initial velocity dispersion) the PS de- 
cays at small scales (large k) as a negative power law in fc, 
while in Hot Dark Matter (hereafter HDM) models (for 
which the prototype is a Universe dominated by a light 
neutrino) there is a exponential cut-off in the spectrum 
(due essentially to the fact that the "hot" neutrinos wipe 
out structures at these scales with their large velocity 
dispersion). All of these models, however, have the same 
"primordial" H-Z spectrum P{k) ~ k (or very close to it) 
on large scales (i.e. small k), that is at scales which are 
large compared to the causal scale at the time of matter- 
radiation equality. This latter scale is of course much 
smaller than our present causal horizon (i.e. than the 
part of the Universe we can probe today). This means, 
in particular, that these primordial density correlations 
should be imprinted in the distribution of matter at very 
large scales, and should in principle be detectable in the 
distribution of galaxies at very large scales, inside the 
present horizon. Until now the only probe of fluctuations 
on such scales is through the temperature variations in 
angle of the CMBR, as the angular correlations in tem- 
perature fluctuations are coupled directly to the three 
dimensional density fluctuations. From the COBE mea- 
surements Q the amplitude of the fluctuations inferred 
is ^ 10^^ in the PS at these scales. We will discuss 
elsewhere the practical difficulties involved in measuring 
such a weak signal in the discrete distribution of galax- 
ies. Here we concentrate on identifying the real space 
properties of these theoretical models at large scales. 

Another context in which an understanding of the sta- 
tistical properties in real space of the H-Z PS of the mass 
density field is important is represented by cosmological 
N-body simulations, the aim of which is to calculate the 
formation of structures under gravity in the Universe by 
a direct numerical calculation (see e.g. (|J^). Because 
the time scale of evolution in these simulations is short 
compared to the dynamical time of the system (i.e. a 
particle moves a small distance relative to the size of the 
box representing a large volume of the Universe) the final 
configuration depends strongly on the initial conditions 
(hereafter IC) at all but the smallest scales. Indeed a 
central idea is that from the final distribution - which 
should be closely related to the observed one of galax- 
ies - one should be able to "reconstruct" some important 
features of the IC, which can be related to other obser- 
vations such as those of the CMBR. A key issue for these 
simulations is thus the setting up of these IC, which in- 
volves subtle problems concerning the discretisation of 
the system. The usual approach to this problem is again 
entirely phrased in fc-space, where instead a real space 
approach proves very useful [Q. To avoid any possible 
confusion for those somewhat familiar with these simu- 
lations, we note here at the outset that the description 



of the H-Z model we give in this paper, as lattice-like or 
glass-like, has no direct relation to the use of lattices or 
glasses in setting up IC in current N-body simulations. 
There lattices or glasses are understood to be sufficiently 
"homogeneous" configurations on which to superimpose 
fluctuations of a desired type. The reason for their use 
(instead of a "uniform" Poisson configuration) is purely 
numerical |Q, and it has nothing to do with the intrin- 
sic statistical properties of the systems being modeled. 
Indeed, as we will explain further at the appropriate be- 
low, these methods have been used primarily to simulate 
cosmological models at the smaller scales at which they 
are not at all glass-like. 

Discussions of real space properties of the density fluc- 
tuations encountered in cosmology are puzzlingly sparse 
in the literature on the subject. Peebles briefly notes ( ^ 
- see pg.523) that a very particular characteristic of H-Z 
models is that "on large scales the fluctuations have to 
be anti-correlated to suppress the root mean square mass 
contrast on the scale of the Hubble length" . Indeed, we 
emphasize the fact that these models are characterized 
at large scales by a correlation function ^(r) which has 
a negative power-law tail: detecting it would be the real 
space equivalent of finding the turnover to H-Z behavior 
to scales at which the PS goes as P{k) ^ k. The prefer- 
ence for a /c-space description is probably rooted in the 
fact that the linear dynamics, which are used to describe 
many problems in cosmology, are most naturally treated 
in this space. While it is true of course that this descrip- 
tion in fc-space is complete, this by no means implies 
that the complementary real space view is redundant, as 
is well known in many contexts in physics. One of the 
points of this paper is to show that this complementary 
view of these apparently so familiar models is at the very 
least interesting and useful, in particular in how it facil- 
itates comparisons with familiar physical systems. 

A basic question we try to answer is the following: 
What "kind" of two-point correlation function is the one 
corresponding to the H-Z behaviour in cosmological mod- 
els? We compare it to some different statistical homo- 
geneous and isotropic systems: (i) Poisson-like distribu- 
tions, (ii) systems with a power-law correlation function 
found in critical phenomena ||l0[| and (iii) distributions 
characterized by long-range order (e.g. lattice or glass- 
like) |ll[]. Through this comparison we can classify H-Z 
models in the third category. We introduce the term 
"super-homogeneous" to refer to these kinds of distribu- 
tions, as their primary characteristic is that mass fluctu- 
ations decay at large scales faster than in a completely 
uncorrelated (Poisson) system. For critical systems one 
has instead a decay of the mass variance which is slower 
than Poisson. Formally the definition of this class of 
"super-homogeneous" distributions is given by the con- 
dition that the PS has P{0) — 0, or equivalently in real 
space that the integral of the two point correlation func- 
tion over all space is zero. In the cosmological litera- 
ture the latter property of cosmological models is often 
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noted, but its meaning (as a strong non-local condition on 
a stochastic process) is not appreciated, or worse misun- 
derstood as a trivial condition applying to any correlated 
system. In the textbook of Padmanabhan , for exam- 
ple, it is "proved" on pg. 171 that the integral over all 
space of the correlation function vanishes independently 
of the functional behaviour of ^ (r) . The error is in an im- 
plicit assumption made that the number of particles in a 
large volume in a single realization converges exactly to 
the ensemble average. This is not true because, in gen- 
eral, extensive quantities such as particle number have 
fluctuations which are increasing functions of the volume 
(e.g. Poissonian, for which the integral is not zero). A 
slightly different, but common, kind of misunderstanding 
of the meaning of the vanishing of the integral over the 
correlation function is evidenced in the book by Kolb & 
Turner ||l^. There it is affirmed (after its statement in 
Eq.(9.39)) to be "...just a statement of mass conserva- 
tion: if galaxies are clustered on small scales, then on 
large scale they must be "anti-clustered" to conserve the 
total amount of mass (number of galaxies)" . The source 
of this misconception seems to be a confusion with the so- 
called "integral constraint" in data analysis (e.g. | |l^ , p^ ), 
which imposes such a condition on the estimator of the 
correlation function in a finite sample, due to the fact 
that the (unknown) average number of points in such a 
sample is estimated by the (exactly known) number of 
points in the actual sample. Despite their apparent sim- 
ilarity, these are different conditions: the first (infinite 
volume) integral constraint provides non-trivial physical 
information about the intrinsic probabilistic nature of 
fluctuations, while the second is just an artifact of the 
boundary conditions which holds in a flnite sample in- 
dependently of the nature of the underlying correlations. 
We will discuss this point in a little more detail at the 
appropriate point below. 

The paper is organized as it follows. In the flrst section 
we recall the basic properties of mass distributions (both 
continuous and discrete) described in terms of stationary 
stochastic processes with a well defined (non-zero) aver- 
age density. In this context we introduce the basic statis- 
tical quantities (homogeneity scale, correlation functions, 
real space mass variance, PS etc.) used to describe these 
systems. We discuss in particular the relation between 
the mass variance in spheres and the PS, noting that for 
power law spectra P{k) ^ /c" and n > 1, the small scale 
(i.e. large k) power dominates the real space variance at 
any scale. We explain that this is not a simple mathemat- 
ical pathology but corresponds to a real property of these 
distributions. Indeed for discrete distributions of points 
we note that a theorem has been proved jl^ showing the 
behaviour approached at n = 1 to be the limiting decay 
of the variance, in real space spheres, in any distribution. 
In the subsequent section we discuss the H-Z spectrum, 
recalling the construction which leads to it in cosmology 
and why it is called "scale-invariant" . We note that the 
H-Z criterion, as naively understood, is not one which is 



satisfied exactly by the spectrum P{k) ~ k. In the next 
section we give a classification of all SSP in terms of the 
behaviour of the PS as fc ^ 0. We give the name "super- 
homogeneous" to those which have P{0) = 0, referring to 
their basic characteristic as more homogeneous than the 
Poisson distribution, with a sub-poissonian decay of their 
mass variance. In the following section we give the exam- 
ples of a lattice, and then a "shuffled" lattice, to illustrate 
the properties of distributions of this type, which have 
the strong order of a lattice or glass at large scales. Here 
we discuss also briefly the relation of our description to 
N-body simulations. In the final section we discuss vari- 
ous points in summarizing our findings. In particular we 
clarify the use of the term "scale-invariance" , "fractal" 
and "correlation length" in relation to the H-Z spectrum 
in cosmology. 

II. BASIC STATISTICAL PROPERTIES AND 
CONCEPTS 

Inhomogeneities in cosmology are described using the 
general framework of stationary stochastic processes 
(hereafter SSP). Let us consider in general the descrip- 
tion of a continuous or a discrete mass distribution p(r) 
in terms of such a process. A stochastic process is com- 
pletely characterized by its "probability density func- 
tional" 7^[p(r)] which gives the probability that the result 
of the stochastic process is the density field p{r) (e.g. see 
Gaussian functional distributions [|l^). For a discrete 
mass distribution the space (e.g. infinite three dimen- 
sional space) is divided into sufficiently small cells and 
the stochastic process consists in occupying or not any 
cell with a point-particle, and p(r) can be written in gen- 
eral as: 

oo 

p(r)=^<5(f-rl), (1) 

i=l 

where fi is the position vector of the particle i of the 
distribution . 

The stationarity refers in the present context to spa- 
tial stationarity of the process, and means that the func- 
tional ^^[^(r)] is invariant under spatial translation. This 
property is also called the statistical homogeneity of the 
distribution. We suppose also that the distribution is 
statistically isotropic (invariance oiV[p{r)] under spatial 
rotation), and has a well defined average value pa, that 
is 

{p{r)) = po > , (2) 

where (...) is the ensemble average over all the possi- 
ble realizations of the stochastic process, i.e. the aver- 
age over the functional V[p{r)]. Statistical homogeneity 
and isotropy (hereafter SHI) imply that the Z-point cor- 
relation functions (p(ri)...p(n)) , for any I, depend only 
on the scalar relative distances among the / points [p^ . 
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Moreover we assume that V[p{r)] is ergodic. In order to 
clarify the meaning of ergodicity, let us take a generic ob- 
servable F = F{p{fi), p{f2), ■■■) of the local density p{f). 
Ergodicity means that (F) is equal to the spatial average 
F given by: 



ro),-) (3) 



where the integral is extended to the whole space fl and 
\\n\\ is its (infinite) volume, and where is (almost) 
any realization of the particle distribution "extracted" 
from the functional V[p{r)]. This property is also referred 
as the self-averaging property of the distribution. Note 
that if the average in Eq. (||) is extended only to a finite 
sub-sample V of the whole space i7 , then Eq. (||) is only 
an estimator of (F) in the given sub-sample. 

In a single realization of the mass distribution the ex- 
istence of a well defined average density implies that 



lim 



d^'r p{r) ^ po > (4) 



oo \\C{R;xo)\\ Jc{R.xo) 

where ||C(i?, xq)!! = 47ri?'^/3 is the volume of a sphere 
C{R,xo) of radius R, centered on the arbitrary point xq 
of space 0. When Eq.| is valid one can then define a 
characteristic homogeneity scale as the scale Aq given by 



'C{R;xo) 



d^rp{f) - po 



< po yR> Ao, Vxo 
(5) 



which depends on the nature of the fluctuations of the 
density in spheres. In practice, in characterizing the scale 
at which a system begins to be homogeneous, it is easier 
to use directly some simple two-point statistics. We will 
mention these definitions at the appropriate point below. 

The quantity {p{'ri)p{f2)...p{fi)) is called the com- 
plete Z-point correlation function. In the discrete 
case {p{fl)p{f"2)---pi'n)) dVi,dV2, ■■■,dVi gives the a pri- 
ori probability of finding I particles, in a single realiza- 
tion, placed in the infinitesimal volumes dVi, dV2, dVi 
respectively around ri,r2, ■■■,fi. 

Let us analyze in further detail the auto-correlation 
properties of these systems. Due to the hypothesis of sta- 
tistical homogeneity and isotropy, {p{rl)p{7^)) depends 
only on ri2 = \fl — f^\. Moreover, (p(rl)p(r5)p(r5)) is 
only a function of ri2 = \fi — t^\, r23 — \f2 — and 
''IS = ~ ''si- The reduced two and three-point corre- 
lation functions ^(r) and C(''i2; ''23, ''13) are respectively 
defined by: 



* Because of the arbitrariness of the position of the center 
of the sphere, the average density is a one-point statistical 
property. 



{p(r\)p(f2)) = pI l + |(ri2) 



(p('-l)p(r2)p(r3)) = 1 + i{ri2) + ^(''23)+ 



^(''13) + C(''12,''23,''13) 



(6) 



(7) 



The correlation function ^(r) is one way to measure the 
" persitence of memory" of spatial variations in the mass 
density Note that, as shown more explicitly below, 
in the discrete case the functions ^ and C differ from the 
usual ^ and C used in cosmology by the so-called diagonal 
part. 

In the discrete case of particle distributions it is very 
important to consider observations from a point occupied 
by a particle. In order to characterize statistically these 
observations it is necessary to define another kind of av- 
erage: the conditional average {F)^^. This is defined as 
an ensemble average with the condition that the origin 
of the coordinates is an occupied point [Q. When only 
one realization p{r) extracted from Vlp^r)] is available, 
can be substituted by the spatial average: 



I N 

F(p(fi),p(f2),...)p = ^E^(p(ri 

i=l 



,Pir2 



(8) 



where the sum is restricted to all the points {N — > 00) 
r'i occupied by a particle of the distribution. Again in 
the case in which the average is restricted to the particles 
r'i belonging to a finite sample of volume V of the whole 
space, we can consider Eq.(||) only as an estimator of 

The quantity 



(p(rl)p(r-5 ) . ..p(n))p dVidV2 ...dVi 



(9) 



gives the average probability of finding I particles placed 
in the infinitesimal volumes dVi, dV2, dVi respectively 
around ri,r2,...,?i with the condition that the ori- 
gin of coordinates is an occupied point. We call 
{p{fl) p{f^) . . . p{f^i)) p conditional /-point density. 

Applying the rules of conditional probability one has 
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{P{r))p 



(p(O)p(r)) 



Po 



(10) 



{p{Q)p{r'{)p{r-2)) 
Po 



However, in general, the following convention is as- 
sumed in the definition of the conditional densities: the 
particle at the origin does not observe itself. Therefore 
{p{r))p is defined only for r > 0, and {p{fi) p(f2)) p for 
fi,T2 > 0. Consequently, and this is what is usually done 
in cosmology 14 1, one can redefine the reduced two and 
three-point correlation function ^(r) and C(''i) ''2, ''12) to 
be equal to ^ and C respectively for r, ri, r2 > 0, and equal 
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to zero for r,ri,r2 = 0. This means simply that the di- 
agonal part is removed from f and ^. In the following we 
use this convention. 

Let us consider the paradigm of a stochastic homo- 
geneous point-mass distribution: the Poisson case. For 
such a particle distribution the reduced two-point corre- 
lation function Eq. m) can be written as (see flS]) 



= — (i.e. CM =0) 
Po 



(11) 



Analogously, one can obtain the three point correlation 
function (Eq.Q): 

C(ri,r2,ri2) = ^^^^i4^ (i.e. C(^i, r2, ns) = 0) . 
Po 

(12) 

The two previous relations are direct consequences of the 
fact that there is no correlation between different spatial 
points. That is, the reduced correlation functions ^ and ^ 
have only the diagonal part. The latter is present in the 
reduced correlation functions of any statistically homoge- 
neous discrete distribution of particles with correlations. 

As already mentioned, in the definition of conditional 
densities, we exclude the contribution of the origin of 
coordinates. Consequently, for a Poisson distribution, 
we obtain from Eq. ([To|): 



{Pir))p = Po 



(13) 



5{f{ - ri) 



PO 



In general ||20|,|l5| for a SHI distribution of particles the 
reduced correlation function can be written as 



Po 



(14) 



6{n)S{ri) 
Pi 



+ Ciri,r2,ri2) 



where ^ and C are the non-diagonal parts which are mean- 
ingful only for r > ad ri , r2 > respectively. In general 
£,{r) is a smooth function of r |2^Jl^. Hence we obtain 
from Eq.|l^ (by excluding again the contribution of the 
origin of coordinates): 



(p(f))p = po[l + e(r)] 
{p{ri)p{r2))p^pl {l + an) 

£.{ri2) + C{ri,r2,ri2) . 



(15) 



-e(r2)- 



be the mass (for a discrete distribution the number of 
particles) inside the sphere C{R) of radius R (and then 



volume II C(i?) II 
is defined as 



-R"^). The normalised mass variance 



a\R) 



{M{Rf) - {M{R)f 

mm? 



where 



{M{R)) = / d\{p{r))=p4C{R)\\. 

Jc{Ft) 



(16) 



(17) 



and 



{M{Rf 



d^n / d^r2{pin)p{f2)) . (18) 

C(fl,) JC{FI.) 

Note that there is no condition on the location of the 
center of the sphere, because of the assumed translational 
invariance of V[p{'r)]. 

In the discrete Poisson case, using Eq. ([ll|), we obtain 
that 

1 1 



a\R) 



(19) 



Po\\CiR)\\ {M{R)) 

In general, for a SHI mass density field with correla- 
tions, substituting Eq.(|^) in Eq.(|l6|), we obtain 

1 



a'iR) = 



\C{RW 



d\i / d\2a\r'i-r'2\) 

C(R) Jc(R) 



(20) 



Using Eq. (|TJ) in the discrete case we can write 
1 



a^R) 



1 



\C{R)\\ 



Po\\C{R)\\ 

C(R) 



C(R) 



d\2^{\r-{~f2\) 



(21) 



Note that the sign of the second term of Eq.(|2l|) is not 
uniquely determined. We clarify this point later on. 
Equations (^^ and (|2^) make evident the relation be- 
tween fluctuations in one-point properties (as in this case 
the number of points in a sphere centered on a random 
point in space) and two-point correlations. In general 
similar links can be found between fluctuations in n-point 
properties and n + 1-point correlations. 

Equation (W) is equivalent to the requirement that 



lim (i?) = 

R—*OG 



(22) 



which is therefore a condition satisfied by any SHI dis- 
tribution. An alternative (slightly different) definition to 
that given by (||) for the scale characterizing homogeneity 
is thus the scale at which cr^{R) reaches unity (or some 
other appropriate fiducial value) []. In the cosmological 



A. The mass variance in a sphere 



In this section we consider the amplitude of the mass tNote that such a definition holds for SHI distributions, and 
fluctuations in a generic sphere of radius R with respect not at all for the case of fractal systems || as discussed in 
to the average mass. First let M{R) = /^,^. p{r)d^r our conclusions section below 
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literature on the distribution of matter (galaxies, clusters 
etc.) in the Universe there is no global convention about 
how this scale is defined; in fact it is a scale which is 
almost never discussed in precise terms. The two most 
commonly used quantities used in characterizing the two- 
point properties are (i) the scale |^ro defined by ^(rg) — 1, 
and (ii) the amplitude of the mass variance at a fiducial 
physical scale, taken to be 8h~^Mpc (e.g ||2^). Given 
(or having determined) the dependence on scale of the 
correlation function or mass variance, these can be easily 
related to simple definitions of the homogeneity scale. A 
practical working definition of homogeneity scale appli- 
cable in the analysis of galaxy surveys, and a discussion 
of the current status of this scale is given in_J§4 25 1. 
Let us return to further discussion of Eqs.(|20|) and (21). 



a2(i?) = 4,,(i?) + S(i?), 



(25) 



It is very important for our discussion to note that this 
condition ( p^ which holds for any mass distribution gen- 
erated by a SSP, is very different from the requirement 



dV^(r) 



(23) 



(where f2 is the whole space) which is a much stronger 
special condition which holds for certain distributions - 
those to which below we will ascribe the name "super- 
homogeneous" . 

Note that, in cosmology (e.g. ^4 15|) the following ap- 
proximation is often used 



dVi / dW(|rl-r-2|)« ||C(i?)|| / dV|(r) 

C{R) JC(K) JC(K) 



(24) 



in particular in evaluating the variance through Eq.(20). 
Such an approximation is not always valid, and the con- 
vergence properties of the double integral need to be ex- 
amined carefully to establish it. In particular it does not 
hold when the condition Eq. ( p3| ) is satisfied. This will 
be evident following the analysis we give below, as we will 
discuss that one has, for any distribution (continuous or 
discrete), a large distance behavior cP{R) = -R"" where 
a < d -|- 1 (where d is the space dimension) . Using the 
approximation ( ^ ) one could apparently obtain through 
Eq.(po|) arbitrarily rapidly decaying behaviors with an 
appropriate power-law behaviour in the correlation func- 
tion. 

In the discrete case, to measure cr^ (i?) one has to take 
into account both terms in Eq.(|2l|), not only the second 
one. From Eq.(|2l|) the variance can, in general, be writ- 
ten as the sum of two contributions: 



■'■This scale has unfortunately been commonly referred to in 
the cosmological literature as the "correlation length" [Q . It 
has no relation to the statistical physics use of the same term, 
which is a scale characterizing the rate of decay of fluctua- 
tions, not their amplitude. See [ pl|j2^ for a clear discussion 
of this point. 



where the first term cTp^j represents the intrinsic Poisson 
noise of any stochastic particle distribution ^, and the 
second term S(r) (which, as noted above, does not have 
to be of a determined sign) is the additional contribution 
due to correlations (i.e. to ^(r) ^ 0). 



B. The power spectrum 

The PS -P(fc) is the main statistical tool used to de- 
scribe cosmological models. It is defined as 



P(k) = (\5,{k)\ 



(26) 



where 5p{k) is the Fourier Transform (hereafter FT) of 
the normalized fluctuation field {p{r) — po)/pQ. For a 
spatially stationary mass distribution p(r) it is possible to 
demonstrate that it can be obtained by simply taking the 
FT of the correlation function ^(r) (up to a multiplicative 
constant) |p6[: 



I 



{271^ 



(27) 



Further, given statistical isotropy P{k) = P{k). For a 
continuous mass density field obtained by a SSP, the two 
basic properties of the PS are the following (Khintchine 
theorem |17 ): 

1) P(fc) > Vfc; 

2) P{k) is integrable in the whole space. 

For a discrete particle distribution the first property is 
still valid, while the second is not because of the diagonal 
part of ^(r) (the Dirac delta function in r = 0). Indeed, 
this part gives a positive constant contribution for ev- 
ery k which makes the integral of P{k) divergent. This 
constant contribution is the PS for the uncorrelated Pois- 
son distribution of particles. Consequently, for discrete 
distributions, the property 2) is modified as follows: 

2') The FT of the ^(r) (i.e. i{r) without the diagonal 
part) is integrable in the whole /c-space. 

In d-dimensions the properties 2) and 2') imply that: 



lim k'^P(k) = 
lim k'^P{k) = 

k — '■oo 



(28) 
(29) 



where in the discrete case P(k) is the FT of (_{r) rather 
than of ^(r). 

In three dimensions we have therefore that, in general, 
the PS can diverge as /c — > with only the condition that 



^Note that this term can give a contribution to the variance 
which dominates over that due to the intrinsic correlations. 
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the divergence is slower that k~^. Any standard type 
cosmological model has in this limit the H-Z spectrum, or 
something close to it, and in any case always has -P(O) = 
0, which implies that Eq.(]23|) holds. We will discuss the 
meaning of this condition at length below. 



C. The PS and real space variance 

Let us analyze the relation between the PS and the 
mass-variance in real space. We first discuss continuous 
density fields, and then make some relevant comments 
on the discrete case. We first rewrite Eqs. gen- 
eralizing them to the case in which we calculate the mass 
variance in a topologically more complex volume V of size 
V. To do this one introduces the window function W\;{f) 
defined as 



Wv{r) = 



1 if f e V 



otherwise 
Therefore we can rewrite Eq. (]l7|) as 



(Af(V)) = / Tyv(rl(p(r))d^r. 
Jn 



(30) 



(31) 



and Eq. (|T|) as 



(m2(V)) - / / d\id''r2Wviri)Wvir2) Wi)p{r2)) , 

(32) 

where the integrals are over all space. The normalised 
variance is then given by 



1 

y2 



d'rid:'r2Wv{?i)Wv{r2)S,{ri ~ ^2) , 



n Jn 



On taking the FT one obtains 
1 



(27r)3 



d^kP{k)\Wv{k)Y 



(33) 



(34) 



which is explicitly positive, and Wv{k) is the FT of 
W\i{r)^ normalised by the volume defined by the window 
function itself, 



Wv{k) = 1 / d^re-'^^-^Wvir) 
^ Jn 



(35) 



with V = J^Wvir)d^r. 

Consider now again the real sphere of radius R for 
which the FT of the window function (normalised as de- 
fined) is 



Wnik) 



{kRY 



(sin kR — kR cos kR) 



(36) 



One then has, assuming statistical isotropy so that 
P{k) = P{k), an expression for the variance in real 
spheres which is 

'^^iR) = 7r^ / dk—--r (sinkR- kRcoskR? k^P(k) . 
27r2 7o (fci?)6 

(37) 

We now show that, for power-law spectra P{k) ~ k" 
(for small k, n > —3) the integral in ( |37| ) has a very dif- 
ferent behaviour for n < 1 and n > 1. For n < 1 the 
integral is dominated by A; ^ R^^, while for n > 1 it be- 
comes dominated by the large k behaviour, and therefore 
sensitive to the PS at a scale k unrelated (in general) to 
R~^. Correspondingly a'^iR) is found to have a limiting 
rapidity of decay at 1/R*, related to the appearance of 
this divergence. Then we give the physical interpreta- 
tion of this result, and note the danger of the use of a 
"Gaussian window" to mask it in cosmology. 

So let us return to Eq.(37) and take a PS P{k) = 
Ak'^e'"^/^" (where A and kc are two constants). We con- 
sider n > — 3 and take the cut-off to satisfy the conver- 
gence properties of the Khintchine theorem. It is easy 
to check subsequently that the results we derive are not 
sensitive to the form of this cut-off at large k. It is con- 
venient to rescale variables to rewrite Eq.(p7|) as 



a^R) = 



1 

27r2 i?3+» 



da: (sin X — a; cos x)^ a;" ''e ^^^^ 

(38) 



with Xc — kcR. 

Since the window function goes to unity when x ^ 
the integral will be well behaved at its lower limit (since 
n > —3). It is also convergent at its upper limit because 
of the exponential. Let us consider the dependence on the 
latter for a sphere much larger than the cut-off scale, i.e. 
R^ k~^, that is Xc ^ 1. For x :s> 1 the integrand goes 
as a;"~2 cos^ x, so that the integral converges, without 
the exponential cut-off, for n < 1. Thus the variance as a 
function of radius behaves as 1/R^^^, and the integral is 
dominated by modes k ^ R^^. In fact, since the integral 
is independent of the cut-off, we have, up to a numerical 
factor of order unity, the relation 



a^R) « -P{k)k%^n- 



(39) 



so that the amplitude of the PS at k can be thought 
simply to correspond to the variance at the physical scale 
R-\ 

For n > 1, on the other-hand, the integral diverges and 
the cut-off comes into play. For n = 1 the integral is 



dx 



In; 



Ini? 



(40) 



so that a'^{R) 
tegral goes as ~ x] 



(In i?)/i?4. Finally for n > 1 the in- 
^ so that one gets a'^{R) ~ l/R"^, 
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independently of n. Importantly, for n > 1, the inte- 
gral in Eq. ( ^8| ) is dominated by the short wavelengths 
with k ~ and not by the fluctuations on the scale 
k ^ R~^, and correspondingly the relation ( |3S| ) does not 
hold. The amplitude of the PS is no longer related to the 
real space fluctuations at the scale k ~ R~^', instead large 
scale spatial fluctuations have their behaviour determined 
by the short scale power in the theory. 

To summarize clearly: For a power-law P{k) ^ fc" 
(with an appropriate cut-off around the wavenumber kc) 
the mass variance for real spheres with radius kc is 
given by 

1. For n < 1, (T^(i?) and the dominant 
contribution comes from the PS modes at A; ^ R^^. 

2. For n > 1, (j'^{R) ^ \/ R^ and the dominant contri- 
bution comes from the PS modes at k~^ . 

3. For n = 1, we have the limiting logarithmic diver- 
gence with a^{R) - (In 

In the cosmological literature j^^] the divergences in the 
latter two cases are treated as a simple mathematical 
pathology due to the assumption of a perfect sphere (with 
a perfectly defined boundary). Replacing the real sphere 
with a smooth Gaussian filter Wv{r) ~ e~'' these 
integrals are also cut-off at the scale k ~ R~^ and one 
recovers a behaviour (J^{R) ~ and a relation of 

the form (|3|) . While of course this is valid mathemati- 
cally it misses an important point, which is that this lim- 
iting behaviour of the variance (as has a very real 
physical meaning which has to do with the nature of sys- 
tems with such a rapidly decaying PS. They correspond 
to extremely homogeneous systems (i.e. extremely or- 
dered systems) in which the variance really is dominated 
by the small scale fluctuations. Let us explain this point 
further. 

Firstly, that the behaviour has nothing to do in prin- 
ciple with the ideality of the perfect sphere is easily seen 
by considering a more realistic modeling of the sphere, 
using a window function giving a smearing on a length 
scale corresponding to the uncertainty in the radius of the 
sphere (This could correspond, for example, to the uncer- 
tainty in the distance measure to a galaxy) . If it is larger 
than the intrinsic cut-off scale in the power spectrum, it 
is this scale which then provides the cut-off in the integral 
giving the mass variance in the sphere. Since this scale 
is in principle independent of the radius of the sphere R, 
the same limiting behaviour of the variance is re- 

covered. Thus it is a physical result for a continuous SSP 
that the mass variance measured in spheres of radius R 
cannot decrease faster than l/i?"*. 



**See, for example, the section entitled "Problems with fil- 
ters" in the book by Lucchin and Coles kA. 



A more intuitive understanding of this fact can be 
gained by considering discrete distributions. One would 
reason that any continuous distribution can be arbitrar- 
ily well approximated at large scales by an appropriate 
discretisation process, and that therefore the same result 
may hold of discrete distributions. In fact such a result 
has been proved several years ago ||2^: In d-dimensions 
there exists no discrete distribution of points in which 
the variance in spheres decays faster than 1/R'^'^^. One 
can see roughly why this is so by considering the most or- 
dered distribution of points one might think of: a simple 
cubic lattice. The variance in a sphere is given by averag- 
ing over spheres with center anywhere in the unit cell. As 
the sphere moves in the unit cell the variance, one would 
guess (correctly!), in the number of points is proportional 
to the difference in the volume of the spheres, which is 
proportional to the surface area of spheres, i.e. oc R'^^^ 
in d-dimensions. Thus the normalised variance scales as 
a result proved in rigorously in |l6| (see also 
[p8| for a more general discussion of the problem) . As we 
will discuss further below the regular lattice, or rather 
a randomized version of it, can be thought of as a kind 
of prototype for the class of distributions to which the 
H-Z spectrum belongs. They are distributions which are 
highly ordered ( "glass-like" ) in which the fluctuations in 
real space actually are at small scales (those at which the 
PS is cut-off). Because of this it is one of their charac- 
teristics, as we have seen, that there is no direct relation 
between the PS at scale k and the physical variance in 
real space at the scale R ^ k^^. 

The Gaussian sphere completely obscures this be- 
haviour for n > 1, giving an apparent behaviour of a 
real space variance oc It does this because it 

models the edge of the sphere as smeared on the length 
scale of the radius (i.e. assumes that the uncertainty in 
our measure of distance to a point is necessarily of or- 
der the distance). Instead of the dependence of the mass 
variance in spheres (with some intrinsic uncertainty in 
the definition of their edges) on the radius, the Gaussian 
sphere gives us the behaviour of the variance in spheres 
as a function of radius when both the radius and the 
smearing imposed on the edge change together and in 
linear proportion. This is a real space measure which 
one can define to recover fc-space properties, but it loses 
completely the essential characteristic of the real space 
variance (measured in real physical spheres). 



8 



III. THE H-Z SPECTRUM AND ITS REAL 
SPACE PROPERTIES 

Let us first recall the kind of argument]^ that sin- 
gles out the H-Z spectrum in cosmology, and why the 
term "scale- invariant" is applied to it. In a homogeneous 
Friedman-Robertson- Walker (hereafter FRW) cosmology 
there is a fundamental characteristic length scale, the 
horizon scale Rnit). It is simply the distance light can 
travel from the Big Bang singularity t — until any given 
time in the evolution of the Universe, and it grows lin- 
early with time. The H-Z criterion can be written 

(tIj{R = Rnit)) = constant, (41) 

i.e it requires that the mass variance at the horizon scale 
be constant. Equivalently, given the proportionality of 
gravitational potential to mass, it can be stated as the 
constraint that the variance in the gravitational poten- 
tial be constant at the horizon scale. It arises naturally 
in the framework of FRW cosmology as a kind of consis- 
tency constraint: the FRW is a cosmological solution for 
a homogeneous Universe, about which fluctuations rep- 
resent an inhomogeneous perturbation. If we take any 
other prescription other than Eq. (^) such a descrip- 
tion will always break down in the past or future, as the 
amplitude of the perturbations become arbitrarily large 
or small. It is in this specific sense that the resulting PS 
is said to be "scale-invariant": there is no characteris- 
tic scale at which fluctuations become large (or small), 
or put another way, they have the same amplitude as a 
function of the only scale in the model. As we will discuss 
further below, it has nothing to do with the same term as 
understood in statistical physics. There scale invariance 
is a characterization not of the amplitude of fluctuations, 
but rather is associated to a particular range of power- 
law behaviors in the correlation function. 

More precisely the form of the H-Z spectrum is arrived 
at from the condition ( p| ) in the following way. We move 
necessarily to a fc space description, as we need to include 
the dynamical evolution of the density field to infer the 
PS inside the horizon today. Let 6k{t) be the amplitude 
of the Fourier component of the density contrast as a 
function of time. To every such mode k we can associate 
a time tc at which it "enters the horizon" , i.e. at which 
the wavelength is equal in size to the horizon. Here 
we work (as almost always in cosmology) with a k which 
is the FT with respect to the spatial coordinates which do 
not change with the expansion, the so-called "comoving" 
coordinates. In these coordinates the time at which the 



mode enters the horizon is given by fcry = 1 where rj is 
the so-called "conformal" time given by ry = / dt/a{t), 
with a{t) the scale factor describing the expansion of all 
physical scales in the Universe. (The horizon scale is 
simply Rnit) = a(t)ri, corresponding to horizon crossing 
criterion {k/a)RH{t) = 1.) The PS today (at t = to, say) 
given by |(5fc(to)P can be written in term of the amplitude 
of each mode k when it entered the horizon. In linear 
perturbation theory, in the matter dominated Universe 
(i.e. recent epochs), the mode evolves as 

4(0=(^)'5.(t). (42) 

In the matter dominated FRW cosmology we have a cx 
t^^^ and thus rj cx t-^^^, so that the time tc{k) when the 
mode k crosses the horizon follows tc{k) cx 1/k^ and 
therefore 

4 (to) cc k^dkitc) . (43) 

The H-Z choice for the primordial PS |(5fc(to)P cx fc is 
then singled out by imposing the criterion 

fc^|4(tc)P = constant , (44) 

which is identified as the mass variance at the horizon 
scale 1] = k~^. We note immediately, following the pre- 
ceding discussion, that the latter identification is in fact 
valid only for power spectra fc" with n < 1. Strictly 
speaking therefore it is impossible to satisfy the H-Z cri- 
terion as it is understood naively; or, to put it another 
way, the H-Z spectrum, that which satisfies (^4|), does 
not satisfy the condition of "scale invariance" since the 
mass variance at the horizon scale (cx rj) is dominated 
in this case by the power at the cut-off scale, not by the 
modes fc ~ rj~^. Taking a spectrum k^~'^ (e > 0) one 
can get arbitrarily close to satisfying the H-Z criterion, 
but the condition of "scale invariance" (in the sense just 
explained) is not physically satisfiable. To avoid this con- 
clusion the criterion could be refined to be that the mass 
variance in Gaussian spheres of radius of the horizon size 
be constant. While it does allow a mathematically co- 
herent formulation, from a physical point of view it is 
an artificial way of avoiding the problem, which is that 
the variance at a given real space scale has nothing to do 
in principle with the amplitude of the PS at the inverse 
scale for n > 1. This is, as we have discussed in the pre- 
vious section, a real physical property of such systems, 
not a mathematical artifact. 

The H-Z spectrum can equivalently be characterised in 
term of fluctuations in the gravitational potential, 64>{f), 
which are linked to the density fluctuations Sp^r) via the 
gravitational Poisson equation PI: 



"We choose here a particular (but commonly used) way of 
describing the HZ spectrum which allows us to avoid too much 
extra formalism. For a commonly used formulation preferred 
by many cosmologists, in terms of a constant "gauge indepen- 
dent" potential, see for example Eg). 



"We simplify here to Newtonian gravity, which becomes a 
good approximation on sub-horizon scales. The comments 
given below can however be generalized to a rigourous formu- 
lation of perturbations to a FRW model. 
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-4:nGSp{r) 



(45) 



From this, transformed to Fourier space, it follows that 
the PS of the potential P^{k) = ^|(50(fc)p\ is related to 
the density PS P{k) as: 

P(fc) 



The H-Z spectrum corresponds therefore to P(j,{k) oc k~'^] 
or, considering the variance in real space spheres of the 
gravitational potential fluctuations, which as for the den- 
sity fluctuations is related to the PS by Eg. , one flnds 
that this variance is constant as a function of R. This 
is the alternative form in which the H-Z condition is of- 
ten formulated. Note that the Khintchinc theorem (cf. 
Eq. (||)) requires that a well defined SSP have P{k) ~ fc" 
with a > —3 for fc ^ 0, so that the H-Z corresponds 
to the limiting (disallowed) behaviour. Equivalently the 
constancy of the variance is in contradiction with Eq. ( |2^ ) 
which requires that the asymptotic variance be zero (in 
order to have a well defined mean about which fiuctua- 
tions are defined) . The H-Z spectrum can thus be seen as 
the (disallowed) limiting behaviour for the potential fluc- 
tuations to be treatable as an SSP. That such a treatment 
be applicable to the potential fluctuations is however not 
a physical requirement. The work of Chandrasekhar H] 
(and see also |^) treats the gravitational force probabil- 
ity distribution in a point distribution and, in particular, 
shows it is well defined even in the Poissonian case, for 
which the potential fluctuations are not an SSP {n = 0). 
To treat the force field as an SSP requires only the weaker 
condition P{k) ~ fc" with a > — 1. 



A. The real space correlation function of 
CDM/HDM models 

All of the current "viable" standard type cosmolog- 
ical models have a "primordial" PS which is the H-Z 
one (or very close to it) down to some arbitrarily small 
scale. During cosmological evolution causal physics mod- 
ifies this spectrum at large fc, which is roughly the causal 
horizon at that time. Around the time at which the mat- 
ter in the Universe (with density scaling as begins 
to dominate over the radiation (with density scaling as 
1/a'*), the evolution becomes purely gravitational at all 
but the very smallest scales, while prior to this time it 
depends strongly on the details of the particular model. 
As a result all such models are H-Z for k < k^q, but 
"turn-over" at this scale to a PS decreasing as a func- 
tion of k. The form of the spectrum in this region de- 
pends on the details of the particular model. Since the 
scale k~^, being the size of the causal horizon at this 
time of matter-radiation equality, is much smaller than 
the causal horizon today, the primordial H-Z PS is in 
principle detectable today. Indirect evidence for its real- 
ity come from the measurements of temperature fiuctua- 
tions in the CMBR, which show a dependence on angular 



scale quite consistent with the H-Z spectrum, with the 
power-law spectrum P{k) fc" giving a fit in the range 
n = 1.1±0.5 Hj. The search to observe this "turn-over" 
to H-Z behaviour directly in three dimensions in the dis- 
tribution of matter at large scales - a central prediction 
and check on such models - has so far proved elusive, be- 
cause of weak statistics at large scales in observations of 
the distribution of galaxies. It is anticipated that forth- 
coming surveys, now being made |Bl| or close to comple- 



tion 1 32 1, will have the capacity to detect this turn-over 
(see the discussion at the conclusions). In the cosmo- 
logical literature this question is again treated almost 
exclusively in k space. Here we look at the character- 
istic real space features which should be found in these 
galaxy surveys if the underlying behaviour is H-Z. In fur- 
ther works we will discuss in detail the question of the 
detection of these features, here we concentrate solely on 
their identification. 

We consider first the two point correlation function. In 
general the FT of the PS of standard cosmological mod- 
els must be done numerically. Before doing so for some 
standard models we considered a simple PS which can be 
transformed, a H-Z spectrum with a simple exponential 
cut-off: 



P{k) ^ Axkx 



(46) 



where A is the amplitude and k^ ^ the cut-off scale. The 
correlation function is given by (see e.g. [|l|) 



A 



For r < Tc = fcf, ^ we have ^(r) 



(47) 



,4 



fc^ > 0, chang- 



ing at r ~ Tc to an asymptotic behaviour ^(r) ~ —r~ . 
Note that the correlation does not oscillate, its only zero 
crossing being at scale r = VSvc- Simply because of the 
condition P{Q) — 0, which implies that the integral of the 
correlation function must be zero, the correlation func- 
tion must change sign and in this case it only does so 
once and thus remains negative at large scales. 

In the normalised mass variance <t^{R) shows a corre- 
sponding change in behaviour from being approximately 
constant at small scales i? < Tc to a In decay at 

large scales, as was shown in Section([IC) above. Note 
that, unlike for the variance in spheres discussed in Sec- 
tion( p^I C| ), there is no limit to the rapidity of the decay of 
the correlation function (for the more general expression 
see iQ ) . Despite the weakness of this correlation at large 
scales, however, the variance in spheres does not behave 
like that of a Poisson system, because of the balance be- 
tween positive correlations at small and negative at large 
scales imposed by the non-local condition P{0) = 0. 

In cosmological HDM models the form of the PS is 
almost the same as we have just considered with an ex- 
ponential cut-off pa 
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with Q = 1,0.2 respectively. The two reference lines have : 

exponents fc , 10° 10' 10' 10' 

r (Mpc) 



P{k) - kexp{~k/kc) 



3/2 



(48) 



A numerical integration verifies that the correlation func- 
tion is essentially unchanged. 

For CDM models, the class by far favored in the last 
few years, the form of the PS at scales below turn-over 
from H-Z behaviour is considerably more complicated. 
In a linear analysis the PS of CDM matter density field 
decays below the turn-over with a power-law ^ fc~^/^ at 
large k until a smaller scale at which it is cut-off with 
an exponential (in a manner similar to that in the HDM 
model). Numerical studies of these models designed to 
include the non-linear evolution bring further modifica- 
tions, roughly increasing the exponent in the negative 
power law regime. For our analysis we have taken an 
analytical approximation to the final PS given by Eisen- 
stien & Hu |3^], and computed numerically the FT to 
the two-point correlation function. We have also com- 
puted directly the variance in spheres. This form of the 
PS is given in terms of the various cosmological param- 
eters. Here we consider for simplicity the case with the 
small baryon density set to zero {fib — 0), which gives a 
PS without the famous oscillations reportedly detected 
in recent observations of the CMBR |3j,^ . This struc- 
ture is not of primary interest to us here because it can 
modify the correlation function only at small scales (it 
arises from causal physics at early times). In Figsj^^we 
show respectively the behaviour of the PS and of the cor- 
relation function for two quite different values of the total 
matter density of the model = 1, 0.2. Minor differences 
will result in the case that there is a cosmological con- 
stant f^A 7^ || . InFig.| we show the behaviour of the 
unconditional variance, computed in real-space spheres. 
We see again a clear convergence in both models to the 
predicted 1/i?'* behaviour beyond the scale characteriz- 
ing the "turn-over" . 

In conclusion two simple real space characteristics to- 
day in the distribution of matter coming from the pri- 
mordial H-Z PS are a negative non-oscillating power-law 



FIG. 2. Behaviour of absolute value of the real space cor- 
relation function for the two CDM models Q = 1,0.2 and 
h — 0.5. The two reference lines are r^* and r~^. Note that 
at small scale ^(r) > 0, with a zero crossing at a scale depend- 
ing on the location of the peak or "turn-over" in the PS, after 
which it remains negative {(,{r) ~ ^J""*) at larger distances. 
The correlation function has been normalised to be ^{ro) = 1 
for ro — 5Mpc. 
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FIG. 3. Behaviour of the unconditional variance in spheres 
for the two CDM models Q = 1,0.2 and h = 0.5. The two 
vertical lines show the transition to the ^(r) ~ behavior 
for the two models. The r~ behavior is a clear and distinctive 
feature corresponding to the P{k) ~ k behavior. 
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tail in the two point correlation function ^(r) ~ —r~'^, 
and a (In R) / decay in the variance of mass in spheres 
of radius R. These are the pecuhar distinctive feature of 
H-Z type spectra which should possibly be detected in 
real space by the new galaxy catalogs. 

IV. REAL SPACE CLASSIFICATION OF LONG 
RANGE FLUCTUATIONS 

We now return to a the discussion of the nature of cor- 
relations in systems with H-Z like power spectra, with the 
aim of elucidating their properties by comparison with 
systems described in statistical physics. To this end we 
first introduce here a classification of all possible mass 
distributions in terms of the main features of the correla- 
tion function ^(r). Following from the discussion of Sec- 
tion||) concerning the behaviour of mass fluctuations, 
we define three distinct classes (for either the case of dis- 
crete particle distribution and of a continuous density 
field): 

1. If 

/ d'^r l{r) = const. > (49) 
Jn 

we can say that at large scale the system is substan- 
tially Poissonian. Indeed Eq. (^ ) implies that the 
PS goes to a constant non-zero value as k goes to 
zero, and therefore that the large distance behavior 
of the mass fluctuations is 

(Af2(i?)) - (Af (i?))2 ^R-^^ {M{R)) . (50) 

We write here the unnormalized form of the vari- 
ance, as the result that the variance of an extensive 
quantity such as the mass is proportional to the vol- 
ume on which it is measured is the most intuitive 
way of characterizing a Poisson type behaviour. In 
this class is, for example, a system with a finite 
range correlation ^(r) ~ e~'^^^''. Beyond the scale 
Tc (the correlation length - see below for a discus- 
sion about this length) the system is uncorrelated 
and effectively Poissonian. 

2. If 

/ d^r ^{r) = -I-C50 (51) 
Jn 

then we are in a case similar to a system at the crit- 
ical point of a second order phase transition (e.g. 
the liquid-gas critical point). Such systems have 
a positive correlation function which is asymptoti- 
cally a positive power law, with i^(r) ~ 1/r'^ and 
7 < d, corresponding to a PS P{k) k''~'^ as 
fc ^ 0. One then has at large scales the variance 



(M^iR)) - {M{R)f - i?" with d<a<2d, 

(52) 

or {M^{R))-(M{R))^ - (MiR))^ with /3 = a/d > 
1. This means that mass fluctuations are large (al- 
ways overwhelming the Poisson fluctuations) and 
thus they are strongly correlated at all scalesj^. It 
is in this context that the concept of self-similarity 
and scale-invariance has been introduced in statis- 
tical mechanics. These terms refer to the fact that 
in these systems the mass fluctuation field has well 
defined fractal properties |^5[] . 

3. If 

[ d'^r |(r) = (53) 

then, as we have discussed, we have for the be- 
haviour of the mass fluctuations 

(M^iR)) - {M{R))^ ~ R°' with d- 1 <a<d, 

(54) 

i.e. (A/2(i?)) - {M{R)f ~ {M{R))P with /3 = 
a/d < 1, so that the mass fluctuations are always 
asymptotically smaller than in the uncorrelated 
Poisson case. This also corresponds to a strongly 
correlated, long-range ordered, system. We will re- 
fer to them with the term "super-homogeneous" 
to underline this feature that they are more ho- 
mogeneous than a Poisson system. (Indeed, the 
Poisson particle distribution is considered as the 
paradigm of a stochastic homogeneous mass distri- 
bution [^). In the context of statistical mechan- 
ics they can be described as glass-like, as they have 
the properties of glasses, which are highly ordered 
compact systems. That can be said to be typically 
lattice-like, with a long-range ordered packing, but 
without the discrete symmetries of an exact lattice. 
Note again that, since ^(0) > (a Dirac delta func- 
tion in the discrete case) by definition, ^(r) must 
change sign with r at least once. They are systems 
with finely balanced positive and negative correla- 
tion. 

The distinction between |l| and |^ is typical of the statis- 
tical physics of critical phenomena in order to distinguish 
a critical state (case ||) from a non-critical state (case 
|l|). In this context the concept of correlation length is 
central. The correlation length is a measure of the dis- 
tance up to which one has spatial memory of the spatial 
variations in the mass density |l9| ]. There is no unique 



''^For example these properties near the critical point of the 
liquid-gas transition gives place to opalescence phenomena. 
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definition of this length scale, but from a phenomenolog- 
ical point of view it can be defined as the length scale up 
to which the effect of a small local perturbation in the 
system is felt. This is due to the fluctuation-dissipation 
theorem which links the response of the system to a lo- 
cal perturbation and the large scale behavior of the two- 
point correlation function (for the different precise defi- 
nitions of the correlation length see for example [|l^). A 
simple definition is (but see also [2l|p^) 



mi 



In case |l| one can generally define a finite correlation 
length, while in case ^ it will generally diverge. In par- 
ticular in the case ^(r) ^ exp(— r/rc), is indeed then 
the correlation length, while for a positive power-law 
^(r) ~ l/r'*' and 7 < d (case ||) Tcoit —>■ 00. 

Case H is typical of ordered compact systems with small 
correlated perturbations. One can meet this kind of cor- 
relation function for example in the statistical physics 
of liquids, glasses, phonons in lattices. The concept of 
correlation length in this context is less central, and the 
extension of its use to this class of systems is not particu- 
larly useful. Instead it is appropriate to classify the cor- 
relation properties of these systems directly through the 
integral of the correlation function as we have done. It 
is this behaviour of their correlations which distinguishes 
them from the other two cases, just as these cases are 
typically distinguished from one another by the value (fi- 
nite or infinite) of their correlation length. Certainly, as 
we have noted, the use of the term "correlation length" 
in the cosmological literature, which is defined as a 
scale defining the amplitude of the correlation function, 
is in no way related to its use in statistical physics. 

Before continuing with a more detailed discussion of 
the nature of this class of super-homogeneous distribu- 
tions to which standard cosmological models belong, we 
clarify one quite common misunderstanding about them 
in cosmology. 



V. P(0) = AND CONSTRAINTS IN A FINITE 
SAMPLE 

As we have noted in the introduction the physical 
meaning of the constraint ^'(0), equivalent to Eq. (53), 



is often missed in the cosmological literature because of a 
confusion with the so-called "integral constraint" , which 
is another very similar, but actually completed different, 
constraint. Let us clarify this point. 

The "integral constraint" refers in this context to a 
constraint which appears in the estimation of the corre- 
lation function in a finite sample (5, say). It is a con- 
straint which indeed can take the superficially similar 
form to (BSh: 



d^r isir) = 



(56) 



where the subscript indicates that the integral is over the 
finite sample volume, and £,e{''') is the value of the esti- 
mator of the correlation function. This is in general a 
quantity calculable from the sample whose ensemble av- 
erage converges to the real correlation function at any 
finite scale when the boundaries of the sample go to in- 
finity. 

That such a constraint has in principle nothing to do 
with the constraint P(0) = is clear from the fact that 
it is one which holds independently of what kind of dis- 
tribution the sample is taken from. Its origin is simple, 
in the fact that the the mean mass density, relative to 
which fluctuations are estimated, is taken in the estima- 
tor from the sample itself. Therefore, roughly speaking, 
the positive correlations measured relative to this density 
are constrained to be balanced by anti-correlations, giv- 
ing rise to a constraint like Eq. (^6|). More specifically 
the two point correlation function can be written as 



{n{r))p 



1 



(57) 



where {n(r))p is the mean density at distance r from an 
occupied point and (n) the true (unconditional) mean 
density. Integrating this expression over the volume of 
the sample S gives the relation 



{Ns)p - {Ns) - (7^) / d'r 
Is 



(58) 



where (A^s)p is the average number of points in the sam- 
ple volume, with a point at the origin by construction, 
while (Ng) is the average number of points is the same 
volume, but without the condition that there is a point 
at the location of the observer. If one estimates the true 
mean density (i.e. (n) = {Ns)/Vs) in a galaxy catalogue 
sample from the actual density in the sample (i.e. on av- 
erage {Ns)p/Vs), the estimator for the correlation func- 
tion will by construction on average obey the condition 
Eq.(|6|) i.e. 



(59) 



Such a systematic (i.e. ensemble average) offset between 
the estimated and the real correlation function is some- 
times referred to in the cosmological literature on the 
subject as "bias" . It is only in very specific circum- 
stances, with certain estimators, that Eq.(p6|) holds for a 
single sample. For an estimator of the form p5[ 



iEir) = 



r(r) 



1 



(60) 



where r(r) is estimator of {n(r))p, and ns an estimate of 
the mean-density from the sample, Eq.(56) will hold if 



ns = ^ I T{r)d\ 
vs Js 



(61) 
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This is in fact a perfectly good prescription for how to 
estimate the mean density in a finite sample, but one 
that is not used in most estimators, which typically have 
a variance around the average behaviour Eq. (|5^). Es- 
timators *** which do not take the mean density as the 
simple density of points in the sample do not in general 
obey even Eq.(p9|), but will always obey some constraint 
of this type, which cannot be avoided because it is intrin- 
sic to the fact that any real sample contains an occupied 
point at the position of the observer. 

In summary there are necessarily constraints on the 
correlation function ^_e(7') measured in a finite sample, 
which may take a form similar to the condition Eq. ( 5^) 
defining super-homogeneous distributions, but over a fi- 
nite integration volume. These two kinds of constraint 
have a completely different origin and meaning, one de- 
scribing an intrinsic property of the fluctuations in a 
certain class of distributions, the other a property of 
the estimated correlation function of any distribution as 
measured in a finite sample. Their formal resemblance 
however is not completely without meaning and can be 
understood as follows: in a super-homogeneous distri- 
bution the fluctuations between samples are extremely 
suppressed, being smaller than poissonian fluctuations; 
in a finite sample a similar behaviour is artificially im- 
posed since one suppresses fluctuations at the scale of 
the sample by construction. An estimator which imposes 
the relation Eq.(p9|) on the estimated correlation func- 
tion would therefore be expected to make a smaller error 
for the class of super- homogeneous distributions than for 
others. We will return to issues such as this in forthcom- 
ing work. 



VI. SUPER-HOMOGENEOUS DISTRIBUTIONS 

In this section we discuss the properties of long-range 
ordered mass distributions, or super-homogeneous distri- 
butions. We first discuss the simplest example of such 
a situation, represented by a lattice of particles. It has 
many of the relevant properties under discussion. By 
studying its perturbation (the "shuffled lattice" ) it is pos- 
sible to understand the properties of more isotropic dis- 
tributions, both continuous and discrete, which are char- 
acterized again by long-range order. The main feature of 
these distributions, as we have discussed, is that a^{R) 
(the unconditional variance) decays faster than in the 
uncorrelated Poissonian case, i.e. faster than i?^'' (where 
d is the space dimension). We then discuss how a con- 
tinuous field with such correlations can be constructed, 
making it clear that the intuitions about the nature of 
the fluctuations in the shuffled lattice can be extended to 



the continuous case. We mention here that one physical 
model in which such correlations are found ^ is the "one 
component plasma" studied, for example, in |ll|j37[ |. This 
models a Coulomb system of discrete positive charges in 
a continuous negatively charged background. In equilib- 
rium the charges reach an extremely ordered glass-like 
configuration with PS at small k like that of the shuffled 
lattice. 



A. The perfect lattice 

The microscopic density in the case of particles placed 
on the sites of a regular lattice (in any dimension) can 
be simply written as 



(62) 



where R is the generic lattice displacement vector and 
ff is position vector of the lattice site with R — Q with 
respect to the origin of coordinates, i.e. R + ff runs over 
all the lattice sites. For simplicity let us suppose we 
have a cubic lattice of unitary lattice spacing. Then, in 
order to eliminate the dependence on the position of the 
origin of coordinates with respect to the lattice, we can 
define an ensemble of lattices by varying the position of 
the origin with uniform probability in an unitary cell. 
For instance in d = 3, —1/2 < r]x,r]y,r]z < 1/2 where 
x,y,z are the axis coordinates, and then the "ensemble 
average" is (...) = / / j^^^^^d-q^dr^ydr^^.... 
Clearly we have that 



W)) = 1 



(63) 



We want to compute now the two-point correlation func- 
tion {p{ri) p{f2)) ■ Firstly we have that 

p{ri)p{r2) = ^(^" - ^)'^(^" R2-rf). (64) 

Ri,R2 

from which one obtains 

{P{ri)p{f2)) = ^ (^^1 - - ^) (65) 

R 

and hence the two-point correlation function 

^>i,r.) = ^^^^p^-l= (66) 



^<5(fi -f2 -i?) -1 



***For a discussion of estimators used in the cosmological 
literature see e.g. |p6|. 



^^^We thank B. Jancovici for describing these systems to us. 
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Note that ^{ri,r2) = £,{fi — r2) which means that our 
occupation stochastic process (i.e. the ensemble) is sta- 
tionary. However since it is not invariant for generic spa- 
tial rotations we have £,{ri — f2) ^ ^2\) (the lattice 
breaks spatial isotropy). 

In order to evaluate P{k) we need to perform the FT 
of ^(r). In the case of a lattice this gives simply 



P(fc)=.^<S(fc-/J) 



(67) 



where the sum is extended to all the dual lattice vector h 
satisfying the duality condition h ■ R = 27rm, where m is 
any integer, but with the exception of /i = 0. Note that, 
because of this last condition, also in this case P{Q) = 0. 

It has been shown ||l^ that for the simple lattice in 
d dimensions the fluctuations in a ballPI of radius R, 



centered on a randomly chosen point, behaves as 



{MiRf) - {M{R))' 



R 



d-l 



(68) 



compared to the Poisson behaviour {M{R)'^) — 
{M{R))^ ^ R'^. This result can be understood as fol- 
lows: 

• In the Poisson distribution, if we take two randomly 
placed spheres of same radius, the numbers of par- 
ticle contained in them differ by an amount which 
is typically of the order of the square root of the 
average number {'-^ R"^); 

• In the case of a lattice the two numbers differ by an 
amount which corresponds to a Poissonian fluctua- 
tion (i.e. the square root) of the number of particles 
contained in the last shell of the sphere of thickness 
equal to the lattice spacing (which scales as i?"^^^). 

This is due to the strong order of the particles in the 
lattice from a large scale point of view. Thus the lat- 
tice has a behaviour of its mass variance which places 
it in the super-homogeneous category, with the limiting 
decreasing behaviour 



a2(i?)^i?-('^+i) , 



(69) 



i.e. a'^{R) - i?"" for d = 3. 



B. The "shuffled" lattice 



shows more evident resemblances with the cosmological 
H-Z case. The recipe is the following: 1) Consider a cubic 
lattice of particle, as discussed above; 2) take a particle 
of the distribution and draw a lattice-oriented cubic box 
of size I larger than the lattice spacing centered on the 
particle itself; 3) displace the particle to a randomly cho- 
sen point of this box; 4) repeat for each particle of the 
lattice. 

We can write p{r) for a certain realization of this 
stochastic process as 



(70) 



where R and fj have the same meaning as before, and ff^ 
is the vector giving the displacement of the particle in 
the box from the lattice rest position R + fj. By defini- 
tion each component of the vector rfg is a random number 
uniformly distributed in the interval [—1/2, -\-l/2]. There- 
fore, for instance in d = 3, the ensemble average (...) is 
now defined to be: 



1 1 1-1/2 ^^^111-1/2 



(71) 



After some algebra one finally obtains that ^{fi, ^2) = 
^{fi — r2) (i-e. the ensemble is stationary) and in d di- 
mensions for integer / one finds exactly: 



|(r) = 6{r) - n 



fc=i 



if Irfel >/ 



< / 



(72) 



It is very simple to verify that J J £_{r)d'^r — 0, 
which is the condition of super-homogeneity. Note that 
for / +00 ^(r) reduces correctly to the simple delta 
function, i.e. to the Poisson correlation function. In fact 
in that limit we must of course obtain a Poissonian distri- 
bution of particles without correlations. Note that conse- 
quently, we have the following non-commutativity of the 
limits: 



= lim 



hm l(r) = 1 . (73) 



One can also find an exact form of P(fc) by applying 
Eq.(p7|) in d-dimensions: 



In this section we define a super-homogeneous stochas- 
tic distribution of particles obtained from a lattice which 



Note that [ |l^ , p8| for the same quantity in cubic boxes 
of size R one obtains {M{Rf) - {M{R)f ~ R'^+^ . This is a 
typical pathology of the lattice which is not a real stochastic 
particle distribution, having a deterministic discrete transla- 
tion symmetry. This pathology is eliminated in the case we 
consider below of a "shuffled lattice" . 



(27r)° 



in 



2(1 - cos(; ki)) 



i=l 



(74) 



Let us analyze the behavior for small values oi k. At the 
leading order we can write 



P{k) 



{2'kY 12 



(75) 



which implies an isotropic behavior for fc — > even 
though by construction P{k) is not isotropic for a general 
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k. Note that we have the P(0) = behaviour of the H-Z 
spectrum. Finally, the fact that for A: ^ oo the PS tends 
to a positive constant, means that such a distribution is 
Poissonian at small scale (r < I). 

It is now easy to calculate analytically the uncondi- 
tional number variance. In particular the calculation can 
be done exactly in cubic boxes with the same lattice sym- 
metry and, with some simple approximations, in spheres. 
In both cases one obtains (J^{R) ~ R~'^ for large scales, 
as in the case of a lattice, but eliminating the pathol- 
ogy of different scaling behaviors between cubic boxes 
and spheres which we noted is present for the rigid lat- 
tice. Note that the fact that (J^{R) ~ at large scales 
corresponds to the fact that, despite the "shuffling" of 
particles with respect to the lattice, the strong lattice or- 
der is maintained at large scales. For non-integer ^ > 1 , 
even though calculations are cumbersome and ^ (r) is not 
simply writtable, the main results about P(k) for small 
k and (7^{R) for large R are the same. 

The H-Z spectrum has this same behaviour character- 
istic of lattice-like order at large scales, while its small 
k PS is P{k) ^ k instead of '--^ fc^. This spectrum cor- 
responds to more power at large scales. We will see in 
the next section that this can be associated with an ap- 
propriately more ordered (i.e. coherent) shuffling of the 
lattice, and precisely what kind of large scale correlations 
is required to obtain the H-Z spectrum will be made ex- 
plicit. The crucial point is that such shuffling must leave 
intact at very large scales the strong order of the lat- 
tice, so that one still has the characteristic behaviors we 
have seen in the shuffled lattice (a correlation function 
which is negative at large scales and integrates to zero, 
a normalised variance in spheres decreasing faster than 
the volume). We thus say that the distribution described 
by the H-Z spectrum has a lattice-like or, more appropri- 
ately because of the isotropy, glass-like long range order. 
More specifically it can be characterized as a glass with 
superimposed opportune coherent long-range perturba- 
tive waves of displacement. 

In relation to this description it is interesting to make 
some brief comments on cosmological N-body simulations 
(see e.g. |^). In this context a standard algorithm used to 
generate initial conditions for these simulations involves 
imposing perturbations on a perfect lattice (or sometimes 
even "glassy" configuration). At first sight this would 
suggest that the point we are making about the H-Z spec- 
trum is in fact already understood in the cosmological 
literature, or at least in the part of it on N-body simula- 
tions. This is not the case, and it is worth explaining this 
to avoid any possible confusion on this point. This tech- 
nique for generating initial conditions has in fact been 
introduced to avoid problems with Poisson noise at small 
scales (large k) in the discretization procedure, and not 
because the system being simulated is understood to ac- 
tually intrinsically resemble a lattice at small k. Indeed 
the primary goal of most of these simulations has been 
to study the dynamical evolution in a range of scales well 



below the "turn-over" in the PS, where the PS (at large 
k) has a negative power-law form (~ k~^, with /3 > 0, 
corresponding to a positive correlation function with a 
"critical" power law behaviour). Thus this procedure is 
applied primarily in a range of scales where the system 
being modeled does not intrinsically resemble a lattice or 
glass at all. 

Only the more recent very large simulations describe 
the larger scales at which the initial conditions should 
have P{k) ~ k. In this case too the use of a displaced 
lattice in setting up initial conditions is not because the 
underlying system is understood to be lattice-like, but is 
simply inherited as a numerical technique for the same 
small scale considerations 30 1. Indeed, as will be further 
clarified in the following section, there is in principal no 
reason why one has to start from a lattice to produce such 
a spectrum; nor indeed is it certain that one obtains the 
right correlation properties if one starts from a lattice. 
What is true is that the spectrum of the initial conditions, 
if it is H-Z, should be glass-like in the sense we have 
discussed. A real space analysis of the initial conditions 
actually used in such simulations shows Q that they do 
not in fact have the appropriate properties. 

Note finally that glass-like systems belong to a wide 
family of distributions for which the common feature is 
that P{k) - k°- with a > for fc ^ and hence P(0) = 0. 
However such behaviors in the PS do not imply directly 
that ^(r) has a negative power-law tail at large scales. In 
particular this is not true if the PS has a singularity for 
P(0) 7^ 0, as happens in many systems. For example Q 
the glass-like distributions (unperturbed and perturbed) 
used as initial conditions in cosmological N-body simu- 
lations have indeed an oscillating ^(r) at all scales, and 



a mass variance cr (r) 



Thus we emphasize that 



the negative power-law tail of the real-space correlation 
function of the H-Z distributions in cosmology is a very 
particular feature of these models. 



C. Uniform distributions with a displacement field 

Let us consider the case of a mass distribution (a den- 
sity field) obtained by superimposing a random displace- 
ment field on a completely uniform density field. 

Let the uniform density field be poir) — po and super- 
impose on it the stochastic displacement field u{r) (the 
infinitesimal volume dV at r is displaced by u{r))- Let us 
call p{r) the resulting density field. We suppose that the 
stochastic field displacement is the realization of a sta- 
tionary and isotropic stochastic process characterized by 
the probability density functional V[u{'r)]. In this way 
V[u{r)] defines also an ensemble of density fields p{r) 
which is stationary and isotropic, with (...) the ensemble 
average. 

By applying the mass conservation (i.e. the continuity 
equation) we find 
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P{r) - Pa 
Po 



(76) 



If we call as usual ^ (r) the reduced two-point correlation 
function of the density field we can write 



C(r) = (V • u(f)V • u{d) 



(77) 



Then, taking the FT of both sides of Eq. (|77|), and mak- 
ing use of the statistical isotropy, we obtain 



P(fc)^fc2p„(fc), 



(78) 



where P{k) is the usual PS of the mass density field and 
Pu{k) is the PS of the displacement field. 

Since Pu{k) is a itself the PS of a SSP it is subject to 
the constraints of the Khintchine theorem. Thus at small 
k it must diverge slower than ^ k~^, allowing to obtain 
at most P{k) ~ k~^, corresponding to a real space corre- 
lation function which must go to zero faster than 1 /r^ as 
r oo. Therefore any continuous SSP of the "substan- 
tially Poissonian" and "super-homogeneous" type can be 
obtained in this way, but not all the "critical" type behav- 
iors. In particular one can obtain a H-Z type spectrum 
with Pu{k) ^ k~^ describing a critical type SSP. 

What is the relation to the discrete case? If we suppose 
that P„(0) — const. > 0, i.e. the displacement field at 
large spatial scales is Poissonian (i.e. uncorrelated), we 
find that for /c ^ one has P{k) ~ fc^. This is exactly 
the same asymptotic behaviour as that we found for the 
case of the shuffled lattice. Indeed we obtained the latter 
through the superposition of an uncorrelated random dis- 
placement field to a "uniform" background, and thus the 
result is natural. Moreover in general we would expect 
the relation Eq. ( |7^ ) to give us in the discrete case the 
large scale behaviour of a set of fluctuations imposed on 
a discretization of the continuous uniform density field, 
and in particular of a rigid lattice which is simply such 
an object. Thus if instead of shuffling the lattice as in the 
previous section we superimpose correlated fluctuations 
with a spectrum ~ k~^ we will obtain at large scales a 
distribution with H-Z behaviour, and the associated real 
space properties at large scales. 

Let us return again flnally to our comments at the end 
of the last section on N-body simulations (sec ||^ for a 
more complete discussion of this point). With the pre- 
vious construction it is easy to see why one arrives at 
the idea of generating a distribution with a certain PS 
by displacing points on a rigid lattice. The lattice is 
simply a discretization of the continuous uniform back- 
ground which is then perturbed. One could in principle 
start from a Poissonian distribution of particles - for the 
generation of the large k behavior - which can be always 
made "uniform" to any desired precision (assuming no 
practical limitation on the number of points used), and 
then superimpose the displacement field to produce the 
required PS. If the distribution produced is to be H-Z, 
and particluary of CDM type, it will have the "super- 
homogeneity" of a lattice at the corresponding scales. 



with the characteristic real space behaviors we have used 
to define it. Equally if one starts from a lattice one can 
arrive at distributions which are not super-homogeneous. 
Indeed, as we noted above, in the context of N-body sim- 
ulations in cosmology the displacement from a lattice to 
produce initial conditions has been introduced in simu- 
lations describing the evolution in cosmological models 
at scales where the models are not H-Z, but rather have 
a positive correlation function with a "critical" power- 
law. The starting point of a lattice has been favored over 
a "uniform" Poisson distribution simply because of nu- 
merical limitations, the latter producing at any feasible 
resolution too much small scale noise overwhelming cor- 
relations at small scales. In summary the central point 
we are making here is not that a H-Z type spectrum can 
be obtained in principal by perturbing a lattice; rather 
the crucial point is that such a system is intrinsically 
lattice-like, irrespective of how a discrete realization of 
it is constructed in practice. These are two completely 
different things. 



VII. DISCUSSION AND CONCLUSIONS 

First we return to the use of several terms "scale- 
invariance" in cosmology. We have described in sec- 



tion([H) with what meaning this term has been intro- 
duced in cosmology: it refers to the fact that the variance 
of the mass (or equivalently gravitational potential) has 
an amplitude at the horizon scale which does not depend 
on time. The PS associated with this behaviour is that 
of a correlated system which is of the super- homogeneous 
type. This use of the term "scale-invariance" therefore is 
not in any way analogous to its (original) use in statis- 
tical physics. In this context it is associated with a dis- 
tinctly different class of distributions which have special 
properties with respect to scale transformations: typi- 
cally critical systems, like a liquid-gas coexistence phase 
at the critical point, which have a well defined homogene- 
ity scale and a reduced two-point correlation function 
which decays as a non-integrable power law: ^(r) r~'^ 
with < 7 < 3. In particular the term does not have 
anything to do with the amplitudes of fluctuations being 
independent of scale: the amplitudes of fluctuations vary 
with scale, while the system is correlated at all scales. 

We note that one might be tempted to associate the 
term "scale-invariance" in the context of cosmology sim- 
ply to the power-law behaviour of the correlation func- 
tion. From a mathematical point of view, the termi- 
nology "scale-invariance" could be used for any distri- 
bution with a power-law tailed £,{r), that is satisfying 
£,{r') = A{b)^{r) for r' — br . In this case the H-Z spec- 
trum would, however, be no more "scale-invariant" than 
any other spectrum with a power-law form at small k. 
In terms of its meaning in physics however this usage is 
restricted to the context of critical phenomena, in which 
it has been demonstrated to be very powerful and use- 
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ful. As we have discussed, this case is completely different 
from the systems we have termed super-homogeneous. In 
particular in the former systems fluctuations are always 
large at all scales, which formally is associated with the 
non-integrability of the correlation function. The super- 
homogeneous systems have correlations of a completely 
different kind: they are delicately balanced to make the 
mass fluctuations smaller than for a Poisson type distri- 
bution. 

In section (IV) we have briefly discussed the use of the 
term "correlation length" . Historically this term was in- 
troduced in the context of critical phenomena to charac- 
terize the transition towards a state of the system where 
(positive) correlations are long range and the normalized 
fluctuations of the mass (or of other extensive quantities) 
decay with the scale more slowly than in the non-critical 
Poissonian state. Moreover in this context, the mean- 
ing of correlation length is given through the fluctuation- 
dissipation theorem, in which the correlation length plays 
the role of the distance up to which the system responds 
to a local external perturbation. It is thus a scale used 
to capture the essential physical distinction between two 
types of distributions which are both in turn distinct from 
the superhomogeneous systems. The essential difference 
between each of these three cases can be best character- 
ized according to the large scale behavior of the mass 
fluctuations in each case. It is therefore useful to take as 
classification parameter the value of the integral of ^(r) 
over all space as we have done, while the notion of corre- 
lation length has no obvious or unique generalization in 
this case. 

Another term whose meaning it is useful to clarify is 
"fractal'. Fractal distributions, which are the prototype 
of scale-invariant geometrical distributions, represent a 
more extreme case of correlated systems: their average 
density is zero in the infinite volume limit and the condi- 
tional average density (^(r))^ decays to zero as a power 
law r~'^ as a function of the distance from an occupied 
point (with < 7 < d). The fractal dimension of the 
mass distribution is given by = d — 7 p^j39| ]. A 
fractal is inhomogeneous at all scales and the concept of 
average density in a finite sample centered on an occu- 
pied point has no intrinsic meaning, because it depends 
on the sample size. Moreover, since the asymptotic av- 
erage density of a fractal distribution is zero, both ^(r) 
and P(k) are undefined for such a system Instead 
one has to work directly with the unnormalized condi- 
tional density {n{r))^. In general, before introducing the 
estimators of ^(r) and P{k) for a finite sample of a sys- 
tem whose underlying properties one does not know (e.g. 
the distribution of galaxies in the Universe), one needs 
to verify that the estimator of the average density po is 
not strongly dependent on the sample size p^ ]. 

There is sometimes confusion in the cosmological liter- 
ature about the meaning of "fractal" in connection with 
the notion of scale invariance. For example in the book by 
Peacock [|o| (Sect. 16) the author writes that "The Zel- 



dovich spectrum is a scale- invariant spectrum....". It is 
then shown that the PS of the fluctuations in the gravita- 
tional potential ^ k^^ and hence the auto-correlation 
function, or the quantity ~ P<t>{k) • fc'^ is a constant: 
the author concludes that "Since potential perturbations 
govern the flatness of space-time, this says that the scale- 
invariant spectrum corresponds to a metric that is frac- 
tal ^: space-time has the same degree of 'wrinkliness' 
of each resolution scale." Both the term "fractal" and 
"scale-invariance" are used here in an incorrect and mis- 
leading way with respect to any of the meanings attached 
to them in their (original) context of statistical physics. 
In particular the H-Z spectrum does not have any prop- 
erties which would allow it to be associated with a fractal 
mass distribution. 

We have discussed the criterion which leads to the H- 
Z spectrum. We have pointed out that the usual naive 
formulation of this condition is in fact not satisfied by 
the spectrum P{k) ~ k, and that one must phrase the 
condition in terms of the variance in Gaussian spheres at 
the scale of the horizon. This is because the spectrum 
P{k) ~ fc is singled out by the constancy (as a function 
of time) of k^P{k) at the horizon scale, which cannot be 
taken to be proportional to the variance in spheres of ra- 
dius R ~ k~^ for P{k) ~ fc" and n > 1. We have empha- 
sized that this is not at all an unphysical behaviour due to 
the ideality of a sphere. Rather it actually tells us some- 
thing very fundamental about the nature of these dis- 
tributions: They are distributions which are so ordered 
at large scales that the variance of mass at large scales 
really does come from small scales. The H-Z spectrum 
marks the transition to a pure lattice-like behaviour of 
the normalised variance in spheres cr^(i?) ^ ^/R^, which 
has been shown to be the most rapid possible decay of 
this quantity for any stochastic distribution of points. 

What then does the use, widespread in cosmology, of 
a Gaussian sphere mean for such a distribution? Mathe- 
matically it simply filters out the power up to a certain 
mode k (which then dominates the integral). Physically 
it can be extremely misleading if interpreted incorrectly 
as a characterization of a variance in real space. Con- 
sider the example of the shuffled lattice. It has PS ^ k^ 
at small k (i.e. kl 1, where I is the shuffling scale). 
Using a Gaussian sphere one would infer that the vari- 
ance at large scale goes as 1/R^. Physically we know 
that all the variance comes from small scales in this case, 
and that this behaviour dominated by the lower cut-off 
gives cr^{R) ^ 1/i?^. The behaviour in Gaussian spheres 
comes from the fact that one is smearing the small scale 
behaviour over the scale k~^. As the sphere grows the 
lower cut-off grows too. One would obtain the same 
behaviour by taking a different smearing scale on the 
sphere's boundary, but by making this scale change in 
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proportion to the radius of the sphere. The behaviour 
observed has to do with the very particular way one is 
smearing, and the real space properties of the system are 
actually obscured. The only usefulness of the Gaussian 
sphere is an alternative way of saying that P{k) ~ k"^, i.e. 
as a statement about k space properties, not real space 
ones. It does not describe in any useful way a property of 
the system in real space. In particular the H-Z criterion 
should be understood really as a A: space one, and cau- 
tion should be applied to its formulation as "constancy 
of mass variance at the horizon scale" . 

We have highlighted the fact that all current cosmo- 
logical models will share at large scales the characteristic 
behaviour in real space of the H-Z spectrum. Specifically 
we note primarily the very characteristic lattice-like be- 
haviour of the variance in spheres cr^{R) ^ (up to 
a small correction which is formally logarithmic for the 
case of exact H-Z), as well as the characteristic negative 
(non-oscillating) power-law tail in the two point correla- 
tion function ^(r) ~ — r~^. In this paper we have not 
addressed practical questions concerning the observation 
of such behaviour in cosmological data. In particular 
one would expect such behaviour to be seen in principle, 
if these models are correct, in the distribution of mat- 
ter in the Universe at large scales, and in particular in 
the distribution of galaxies. So far such behaviour has 
not been observed. Rather the characteristic feature of 
galaxy clustering at small scales is that it shows fractal 
behaviour (2^,^l[l , which as we have noted corresponds to 
a very different kind of distribution than that described 
by CDM type models. A centrally important (and much 
debated |2^,^,|2^) observational question is the deter- 
mination of the scale marking the transition from this 
behaviour to homogeneity. On much larger scales galaxy 
structures should then present the super-homogeneous 
character of the H-Z type PS. Indeed this should be 
a critical test of the interpretation of measurements of 
CMBR in terms of the H-Z picture on large spatial scales 
P, ^4|j3^ . Observationally a crucial question is the feasi- 
bility of measuring the transition between these regimes 
directly in galaxy distributions. With large forthcom- 
ing galaxy surveys it may be possible to do so, but this 
is a question which must address exactly the statistics 
of these surveys and the exact nature of the signal in 
any given model. These are questions we will address in 
future works. One other direct usage of the results de- 
veloped here is in the context of N-body simulations, in 
which one studies numerically the evolution of perturba- 
tions in cosmological models, and a knowledge of their 
real space characteristics can be very useful. We refer to 
for a detailed discussion of this point. 
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